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We have studied the superconducting state of Sr2Ru04 under a magnetic field parallel to 
the superconducting plane. We show that due to a weak spin-orbit coupling, non-unitary 
k y (z — iay) state is formed right at H c2 , then changes to unitary k y z state, as a magnetic 
field is lowered. In terms of this crossover, we address the origin of the observed double peaks 
of specific heat and the disappearance of the double peaks at low fields. 
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Sr 2 Ru04 has been intensively studied in the last years, and its superconducting properties 
are well understood as to a zero magnetic field state. 1 ' 2 /uSR experiment 3 and a microscopic 
calculation 4, 5 consistently support the chiral triplet superconductivity with the order param- 
eter, denoted as d(k) oc (k x ± ik y )z. Recently, momentum dependence of the gap amplitude 
was also determined with a high-precision measurement of the angle-resolved specific heat. 6 

However, the two-component order parameter d(k) oc {k x ± ik y )z seems to contradict a 
behavior under the in-plane magnetic field. As Agterberg pointed out, the in-plane magnetic 
field lifts the degeneracy between k x i and k y z. Given a magnetic field is applied parallel to 
the x-axis, only one component with a shorter coherence length, say kyi, is condensed near the 
upper-critical field (H C 2)- Since k y z is symmetric with respect to a reflection about the y-axis 
while (k x ± ik y )z are not, a phase transition line must exist between H = and H = H C 2 for 
all temperatures below T c . 

Actually, the observed double peaks of the specific heat, 8 ' 9 the thermal conductivity, 10 
and the ac-susceptibility 11 suggest an existence of a second-order phase transition. However, 
these observations show several inconsistencies with the prediction by Agterberg. 7 Firstly, 
the second peak was observed at much higher field than the theoretical estimate. Secondly, 
the second peak line terminates at (H, T)~(1.2T, 0.8K) and merges into the line. In 
particular, the termination of the second peak line suggests that the zero-field (k x db ik y )z 
state can be smoothly continued up to H C 2 near T c , thus inconsistent with the symmetry 
argument above. 

Since the position of transition line is susceptible to the form of the pairing function, 
etc., it may be possible to attribute the observed peaks to Agterberg's phase transition. 12 
However, in this paper, we would like to address this problem from a different standpoint, 
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i.e., by considering the d-vector degrees of freedom. In Sr2Ru04, at zero field, d-vector is 
fixed to be parallel to z due to the spin-orbit interaction between 3d electrons forming the 
three Fermi surfaces. 5 However, the energy splittings between the stable (A)(k x ± ik y )i and 
unstable (B)k x x± k y y or (C)k x y±k y x configurations are considerably small. According to the 
microscopic calculation, it was estimated to be less than 0.01T C . 5 Such a small energy difference 
naturally results from the fact that spin-orbit interaction affects the superconductivity only 
in the order of (pp) 2 , where A is the LS coupling constant, and W is the band-width. These 
tiny energy differences are also supported by the NMR experiment with the magnetic field 
parallel to z. 13 There the d-vector is shown to flip at H ~ 200Oe, which means the states 
(A)-(C) are degenerate within ~0.001T C . 

First, let us explain qualitatively the physical mechanism of the double peaks which is 
proposed in this paper. The quasi-degeneracy of order parameters enables the k y y component 
to mix with the most stable k y z. Actually, under a magnetic field, the mixing is caused by 
the spin-dependent chemical potential shift due to Zeeman coupling with the external field. If 
the density of states(DOS) becomes larger with increasing energy as in Sr2Ru04, the Cooper 
pairs with spins parallel to H are easier to form than those with anti-parallel spins. This effect 
is especially strong in Sr2Ru04 due to the Van Hove singularity just above the Fermi energy. 
Therefore, the spin-polarized Cooper pairs, i.e., k y {z — iay) state becomes stable near H C 2, at 
the expense of part of spin-orbit coupling energy. As the magnetic field is lowered, however, 
the k y (z — iay) state becomes unstable, since a non-unitary state generally loses condensation 
energy. As a result, a crossover to the unitary k y z state occurs accompanying a large entropy 
release. This crossover naturally leads to the double peaks of specific heat as observed. 

We also explain the disappearance of double peaks near T = T c at low magnetic fields 
in the same mechanism. The crossover is observable only when the energy gain due to 
spin-polarization of Cooper pairs surpasses the spin-orbit coupling energy. While the spin- 
polarization becomes stronger with increasing magnetic field, the spin-orbit coupling energy 
does not change with field. Therefore, the double peaks are observed only at low temperatures 
when H C 2 is high enough. 

In the following, we demonstrate our ideas given above by using the following effective 
Hamiltonian, 




E [(**(k)^(k') + * y (k)* y (k'))ct )S ci kia c_ lc , ia c kfia 



k,k',s 



\ 3 2 E ["z(^(k)* w (k')-^(k)* x (k'))ci iS C7i k)a C'_ kf)a C' k , )a 
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^(k)^(k') - ^(k)^(k'))4,T^lk,T C -k'a^kM + /i.e. . (1) 

k,k' " 

Here Ck, s represents electron annihilation operator for 7-band. s is the Pauli spin matrix and 
|(|) means the spin direction (anti)parallel to z. The pairing functions 3> x (k) and ^(k), 
and the coupling constants g±-g4 are obtained from the results in the third-order pertur- 
bation calculation on a Hubbard-like model. 5 <fr x (k) and ^(k) have the same rotational 
properties as k x and k y under the point group of Sr2Ru04. The detailed forms of these 
function do not matter in the subsequent analysis. The kinetic energy term £k is expressed 
with the two-dimensional tight-binding model for the 7-band, £k = — 2t z (cos(k x ) + cos(k y )) — 
4t' z cos(k x ) cos (k y ) — [x z , with ^ = 0.4 and (i z ~ 1.555. Note that we are interested in a 
sufficiently high magnetic field where contributions from the a- and (3- bands are negligible. 

Values of gj are determined in order to reproduce the transition temperature, T c , of each 
order parameter. T c of each order parameter component is given by the combination of gj as 
shown in Table. I. When the spin-orbit coupling is absent, i.e., A = 0, all the states (A)-(C) 
are degenerate, i.e., g\ = 53 and 54 = 0. However, the spin-orbit coupling lifts this degeneracy 
and (A)(3> x (k) ± ^^(k))^ shows the highest T c , then 53 becomes dominant. g% — g\ > and 
54 are originated from the spin-orbit interaction and much smaller than g$. 54 determines 
the relative stability of (B)$ x (k)x ± 3> y (k)y and (C)$ x (k)y ± <& y (k)x. Although the two-fold 
degeneracies in (B) and in (C) are lifted by a weak interaction through the #2 channel, #2 
will play no role in the subsequent analysis and we put g% = for simplicity. Below, we use 

NW)Fi ~ NWE = °- 001 and Af(0)(gi+ g4 ) - jV(0)(g!-g4) = °- 0015 as *yP ical Parameters. 





pairing function 


G 


(A) 


($ x (k) ± i$„(k))z 


g3 


(B) 


<S> x (k)k±%(k)y 


gl-g4 


(C) 


$ x (k)y±$ y (k)x 


gl+g4 



Table I. Classification of pairing functions and the corresponding coupling constants. The transition 
temperature of each state is proportional to exp[— j^jq ] with N(e) being the density of states at 
e = ep. 



In order to study the superconducting state close to H C 2, we derive Ginzburg-Landau 
free energy from eq. (1) in the weak-coupling approximation, as is legitimate for Sr2Ru04 
(T c /ep ~ 10~ 4 ). We ignore the d-vector component parallel to H = (H cos 4>o, H sin </>o, 0), 
which is suppressed due to the Pauli-paramagnetic effect. Furthermore, we only consider 
the components with short coherence lengths, i.e., those with the pairing function <J>||(k) = 
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Fig. 1. The direction of magnetic field H. The direction parallel to H is termed as ||. The direction 
perpendicular to both H and z is termed as _L. 



cos0o^i/(k) + sin ^^(k). 14 By taking the spin axis parallel to H, we obtain 



= K\ d ± \ 2 + \ d z | 2 ) + a so | d ± | 2 +C 2± [| d ± d ± | 2 + | d ± d z | 2 ] + C 2z [| 8 z d ± | 2 + | 8 z d z | 2 ] 

2 



ia sp (d±d* z - d z d\) + ^ (| d± | 2 + | d z | 2 ) 2 + | (d±d* z - d z d* ± ) 



1 / 1 1 \ 

( ' >o TT7777" ( — ■ x J, (3) 



iV(0) v « ? i + (74COS 2 (2 ( / ) o) 53 
1 f 1 rr . nr/ . 1 tanh(^) 



/l 1 tanhf — ) 

+ -<^b#) - iV(£ k - 2^fl^)] ^ (> 0), (4) 

where the order parameter dj_ and d z are defined as 

d± = -(A n + A u ), (5) 

4 = ^(A TT -A U ). (6) 

Here, Af-pm) means the order parameter component for the electron pairs with (anti) parallel 
spins to the magnetic field H. The gauge invariant differential operator is given by dj = 
Vj — i-§Aj (j =_L,z), where the subscript _L means the direction perpendicular to both 
H and z (see Fig. 1). The coefficients in eq. (2) are t = T j^° with T c q being the transition 
temperature at H = for (A), and C 2z{±) = ^^(^^(k) 2 )^, C 4 = ig%<$|| (k) 4 )^ 
with (•■ - )fs being the average over the Fermi surface. 

The two parameters a so and a sp play important roles in the subsequent analysis. a so 
is originated from the spin-orbit interaction, which stabilizes d z compared with d±. On the 
other hand, a sp stems from the spin-polarizing effect, which favors Atj compared with An - To 
understand this, it is helpful to rewrite this term as — ia sp (d±d* — d z d* ± ) = — ^a sp (\ Aj| | 2 — | 
Aj| | 2 ). As explained in the Introduction, this term comes from the spin-dependent chemical 
potential shift due to the gradient of density of states. For the spatial dependence of the 
order parameters, we introduce the parameters r, 8, and Xi an d put d^ = rsm9^(r) and 
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dj_ = r cos 9e~ lx ^> (r) (0 < 9 < ^), where \&(r) is the wave function in the lowest Landau 
level, satisfying the relation {C 2 ±_d]_ + C 2z dl)^(r) = VC 2 _lC 2z ^H^(t) = ^*(r) and being 
normalized as J * | ^(r) | 2 = 1. Here, H$ is the upper-critical field at T = when a sp = 
0. Considering that Hq ~ 1.5 Tesla for Sr2Ru04, we set g ^ H ° = 0.5. Using these order 
parameters, we obtain the spatial average of Fql {Fql(t,9,x)) and minimize Fql with respect 
to r and \- The results are 



(t c - t - ^( 2a sp y + al S m 2 (6 - 9 c )f 
l + sin 2 (2#) 



FGL{r m in, 9, Xmin) OC | — j— , (7) 



where 



2 | t | -a so cos 2 6> + a sp sin(26>) - 
Tmm = C 4 | ^ | 4 (l + sin 2 (20)) ' (8) 

Xmin = ~^ ; (9) 



<c=-;t + v a - +( ^ )2 -^' (11) 

The critical temperature t = t c was found by putting r 2 = and 9 = 9 C in eq. (8). Here, # c 
expresses the spin-polarization of Cooper pairs at t = t c , and its value is determined from 
the competition between a so and a sp as in eq. (10). If spin-orbit coupling is strong enough 
(a so » a sp ), 9 C becomes ~ |, i.e., the state approaches a unitary state $||(k)z. Whereas, if 
spin-orbit coupling is sufficiently weak (a so < a sp ), 9 C deviates from ^ and then a non- unitary 
state <fr||(k)(z + iay) is formed. Below t c , free energy is obtained by minimizing eq. (7) with 
respect to 9. From eq. (7), one can see that 9 approaches ^ much below t c , and the unitary 
state <3?||(k)z is formed irrespective of the value of a sp . 

From the minimized free energy, we obtain the temperature dependence of the specific heat 
C for various fields, which is shown in Fig. 2. Here, in order to obtain C correctly to the order of 
t c — t, we have perturbatively included the sixth order term ^^^t^ - | ^ | 6 r 6 (l + 3sin 2 29) 
in the free energy. As shown in Fig. 2, C has two peaks at high magnetic fields, while only a 
single peak at low fields, consistent with the experiments. In order to characterize the crossover 
region, we define a crossover temperature t* and a crossover magnetic field H* where C takes 
the minimum. t*s are shown with arrows in Fig. 2. t*(H*) can be associated with the starting 
point of the crossover with decreasing temperature (field). In Fig. 3, we plot H*(T) in T- 
H diagram. H*(T) line is located just below H c2 , and is terminated at T ~ T t = 0.50T c o. 
Since the double peak structure becomes unobservable below T t , this point can be naturally 
identified with the observed "tricritical point" as shown by an arrow in Fig. 3. This point 
corresponds to (T t , i^*)=(1.2K, 0.8T). 
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As to the field-angle dependence of the crossover behavior, we could not find any marked 
differences between <j)Q = and (f>o = j- Actually, while the ac-susceptibility shows the charac- 
teristic field-angle dependence, 8,9 the double-peak structure of specific heat does not change 
with field direction, thus consistent with our analysis. 

Let us make some additional remarks. Firstly, the crossover discussed in this paper has a 
lot of features in common with the Ai— A phase transition in 3 He under magnetic field. 15-17 
Due to a lack of spin-orbit coupling, a phase transition occurs in 3 He instead of a crossover. In 
order to discuss the Ai— A transition of 3 He in our context, one only has to put a so = in eq. 
(2). Then, one can find that the Ai— A transition takes place at t = t c — 2a sp . Secondly, our 
result does not deny the Agterberg's transition 7 at low magnetic fields. Since the k y (z — iay) 
state cannot be adiabatically continued to (k x ± ik y )z state, there still has to be a second- 
order phase transition at a finite magnetic field up to T=T c o. Actually, it was found that 
magnetization shows a kink at low magnetic field, suggestive of the phase transition. 18 

Thirdly, we would like to comment on the zero-field state. So far, we have assumed that 
the (k x ±ik y )z state appears at zero magnetic field. However, this assumption is not necessary. 
For example, if the k x x±k y y is stable and is nearly degenerate with (k x ±ik y )z, one will find 
a crossover from k y (y + iaz) to k y y, and double peaks appear in specific heat in the same 
way. Finally, the second-peak line of specific heat is placed at H ~ 0.95ff C 2 at T ~ 0, which 
is lower than our theoretical estimate H* ~ 0.99i/ C 2. However, in our analysis, we have not 
considered the renormalization of the effective mass(m*) and g- factor (g*). Due to the electron 
correlation, these factors are enhanced compared with the bare value as — ^ — ~ 5.5 and 
— — ~ 1.5, respectively. Since a sp is proportional to the product m*g*, the renormalization 

Qbare • r 

of these values will push the crossover line down to lower field. Hence, the observed second- 
peak line can be reproduced even quantitatively. 

When a magnetic field is applied slightly off the plane, it is observed that the double 
peaks of the specific heat disappear. 9 Currently, our theory cannot explain this phenomenon. 
However, the deviation of H c2 line from the result of Ginzburg-Landau anisotropic effective 
mass model implies that we need to take account of the layered structure of this material in 
order to reproduce this phenomenon. We would like to leave this problem for future task. 

In summary, we have shown that the observed double peaks of specific heat in Sr2Ru04 
can be attributed to the crossover from k y (z — iay) to k y z state. Within our theory, we could 
also explain the disappearance of double peaks at low magnetic fields. 

We would like to thank M. Ichioka, M. Arai, and S. Yoshio for valuable discussions. 
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Fig. 2. Temperature dependence of the specific heat for =0.20, 0.40, 0.60, 0.80, 0.90 from right 
to left, y is normalized so that the jump at t = t c is 1 when a sp = 0. The crossover temperature, 
t* , are shown with arrows for -f- =0.60, 0.80, and 0.90. For other fields, t* is absent. The inset 
shows the H dependence of 6 C . 




Fig. 3. H*(T) line is shown with a red line in the T-H diagram. A black line shows iJ C 2 line. The 
"tricritical point", {T t , H t ), is shown with an arrow. (Inset) Same with the main panel, with H 
normalized with H C 2 of each temperature. 
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